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We present measurements of the time-dependent CP- violating asymmetries in _B° K^K^Ks 
decays based on 384 million T(4S') BB decays collected with the BABAR detector at the PEP- 
II asymmetric-energy B Factory at SLAC. We obtain the CP asymmetry parameters C = 0.02 ± 
0.21 ± 0.05 and S — —0.71 ± 0.24 ± 0.04, where the first uncertainties are statistical and the second 
systematic. These results are consistent with standard model expectations. 



PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 

In the standard model (SM) of particle physics, the de- 
cays i?" ifgi^gifg are dominated by the b sss glu- 
onic penguin amplitude. A large violation of CP symme- 
try is predicted by the SM in the proper-time dependence 
of b CCS decays of neutral B mesons. Recent measure- 
ments of CP violation in 6 — > ccs decays 1] are in good 
agreement with the SM prediction [4]. The predicted 
amplitude of this CP violation (CPV) is sin 2/3, where 
/? = aj:g{-VcdV*^/VtdVfl) is defined in terms of the ele- 
ments Vij of the Cabibbo-Kobayashi-Maskawa (CKM) Q 
quark mixing matrix. The SM also predicts that the am- 
plitude of time-dependent CPV in b —^ sqq (q = d, s) de- 
cays, defined as sin 2/3cff, is approximately equal to sin 2/3. 
Contributions from loops involving non-SM particles can 
give large corrections to the time-dependent CPV am- 
plitudes for these decays. The theoretical uncertainty in 
the SM prediction of sin2/3off is particularly small, less 
than 4%, for the decay B° K'^K^K^, which is a pure 
CP-even eigenstate A violation of sin2/3eff — sin 2/3 
would be a clear sign of physics beyond the SM Q . In this 
paper we present a measurement of the time-dependent 
CP- violating asymmetries in the decay i?" K'^K';!K'^ 
0. 



The results presented here are based on 383.6 ± 4.2 
milhon 7(45*) -> BB decays collected with the BaBAR 
detector at the PEP-II asymmetric-energy e+e^ collider, 
located at the Stanford Linear Accelerator Center. The 
BaBAR detector [Tl] measures the trajectories of charged 
particles with a five-layer double-sided silicon microstrip 
detector (SVT) and a 40-layer central drift chamber 
(DCH), both operating in a uniform 1.5 T magnetic field. 
Charged kaons and pions are identified using measure- 
ments of particle energy-loss in the SVT and DCH, and 
of the Cherenkov cone angle in a detector of internally 
reflected Cherenkov light. A segmented CsI(Tl) electro- 
magnetic calorimeter (EMC) provides photon detection 
and electron identification. Finally, the instrumented 
flux return of the magnet allows discrimination of muons 
from pions. 



The time-dependent CP asymmetries are functions of 
the proper-time difference At = tcp — itag between a fully 
reconstructed P° K^K^K2; decay {Bcp) and the other 
B meson decay in the event (Stag), which is partially 
reconstructed. The decay rate /+ (/_) when the tagging 



meson is a B^ [B^] is given as 

p-IAtl/r^o 

/±(Ai) « — X (1) 

4rB0 

[1 ± S'sin(AmrfAt) =F C'cos(AmdAt) ] , 

where tb^ is the i?" lifetime and Am^ is the B'^-B^ mix- 
ing frequency. The parameters C and S describe the 
amount of CP violation in decay and in the interference 
between decays with and without mixing, respectively. 
Neglecting CKM-suppressed decay amplitudes, we ex- 
pect S" = - sin 2/3 and C in the SM. 

The data are divided into two subsamples, one where 
all three K'^ mesons decay into the tt+tt" channel 
{Bcp{-\ )) and another where one of the if^ mesons de- 
cays into the tt^tt^ channel, while the other two decay 
into the tt^tt^ channel (i?cP(oo))- 

We form 7r° 77 candidates from pairs of photon 
candidates in the EMC. An energy deposit in the EMC 
is determined to be a photon candidate if no track in- 
tersects any of its crystals, it has a minimum energy of 
50MeV, and it has the expected lateral shower shape 
in the EMC. We reconstruct tt'^tt'^ candidates 

from tt'^ pairs with an invariant mass in the range 480 < 
m^o^o < 520MeV/c^. We reconstruct tt+tt" can- 

didates from pairs of oppositely charged tracks, originat- 
ing from a common vertex, with an invariant mass within 
12 MeV/c^ (about 4 standard deviations) of the nominal 
mass |2[. We also require the decay vertex to be 
along the expected flight path and the significance of the 
reconstructed flight distance t^o /ctt- □ to be larger than 

S Kg 

5. 

For each B(jp(^_^ ) candidate two nearly independent 

kinematic variables are computed; the beam-energy- 
substituted mass Toes — \/ {s/2 + Pj • pbY /E'f — p-^, 
and the energy difference A_E = — Vs/2. Here, 
(-Ei,Pi) = qe+e- is the four-momentum of the initial 
e^e~ system in the laboratory frame and is the 
center-of-mass energy, while ps is the reconstructed mo- 
mentum of the B^ candidate in the laboratory frame and 
E*g is its energy calculated in the e^e~ rest frame. For 
each i?cP(oo) candidate we use two different kinematic 
variables. They are the reconstructed P" mass tub and 
the missing mass TOmiss — \/ {qe+e- — to)^, where is 
the four-momentum of the i3cP(oo) candidate after a mass 
constraint on the P° meson has been applied. Due to 
leakage effects in the EMC, which affect the photon en- 
ergy measurement and therefore the tt*^ reconstruction, 
the shape of the tub distribution is asymmetric around 
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the mean value. This results in this combination of vari- 
ables being less correlated than AE and mEs , with better 
background suppression Q . 

For Bcp signal decays, the mEs, Wmiss and uib dis- 
tributions peak near the i?" mass, while the AE distri- 
bution peaks near zero. For Bcp{-\ ) candidates, we re- 
quire 5.22 < TOES < 5.30 GeV/c2 and |A£;| < 120 MeV. 
For Bcp{oo) candidates, we require 5.11 < m,niss < 
5.31 GeV/c^ and |tob - rn^^^\ < 150 MeV/c^, where 
nig^'^ represents the world-average i?° mass Q . These 
selection windows include the signal peak and a "side- 
band" region which is used for characterization of the 
background. 

The sample of Bcp candidates is dominated by random 
KgKgK'^ combinations from e^e^ ~^ <N {<1 — u, d, s, c) 
fragmentation (the qq continuum). We use topological 
observables to discriminate jet-like e+e~ qq events 
from the more spherical BB events. In the e"'"e~ rest 
frame we compute the angle 0^ between the thrust axis 

of the B(jp^_^ ■) (50^(00)) candidate's decay products and 

that of the remaining particles in the event. We require 
I cos < 0.90(0.95), which reduces the number of back- 
ground events by one order of magnitude. We also use 
the Legendre monomials Ln and L2 , for the characteri- 
zation of the event shape [8|. The monomials are com- 
bined in a Fisher discriminant 3 (ratio I2 = L2/LQ) 

for i3cp(H ) {Bcp(Qo)) candidates, and it is used in the 

maximum-likelihood fit described below. 

The average Bqp candidate multiplicity in the -Bc'P(oo) 
sample is approximately 1.7, coming from multiple 
tt'^tt'^ combinations. In these events, we select the com- 
bination with the smallest = ~ "^Kl)"^ / '^Im^ 
where rrii (m^^o) is the measured (world-average) 
mass [1] and Orm is its estimated uncertainty. We use the 

same method in the Bcp[^ ) sample, where only 1.4% 

of events have more than one B(jp(^^ ■) candidate. 

Since 5° decays proceed through a 6 ^ ccs 

transition, we remove all Bcp^_^ )(Scp(oo)) candidates 

with a K'^K'^ mass combination within 3cr {2a) of the Xco 
or Xc2 mass. After these vetoes, the total reconstruction 
efficiency, including K'^ branching fractions, is about 
6% (3%) for Bcp(+-) (BcP(oo)) candidates, assuming a 
uniform Dalitz distribution. 

The remaining background from BB events is esti- 
mated to be negligible for the B(jp(^^ ) sample and is 

absorbed into the qq continuum component. For the 
^CP(oo) sample, we extract the yield of BB background 
events simultaneously with the signal and qq event yields. 

A multivariate tagging algorithm determines the fla- 
vor of the Stag meson and classifies it in one of seven 
mutually exclusive tagging categories [J Q. They rely 
upon the presence of prompt leptons, or one or more 
charged kaons and pious in the event, and have differ- 
ent purities. We measure the performance of this algo- 
rithm with a data sample (-Bfiav) of fully reconstructed 



B" 7r+/p+/af decays. The effective tagging ef- 

ficiency is Q = ^^£^=(1 - 2w''f = 0.304 ± 0.003, where 
e'^ {w'^) is the efficiency (mistag probability) for events 
tagged in category c. 

We compute the proper-time difference At — Az/j(3c 
using the known boost of the e^e^ system and the mea- 
sured separation between the Bcp and Stag decay ver- 
tices along the boost direction (A2 — zcp — ^tag) 0. 
For the Bcp decay, where no charged particles are pro- 
duced at the decay vertex, we determine the decay point 
by constraining the B production vertex to the interac- 
tion point (IP) in the plane orthogonal to the beam axis 
using only the — > tt+tt" trajectories. The IP posi- 
tion is determined on a run-by-run basis from two-track 
events. We compute At and its uncertainty (TAt from 
a geometric fit to the T(4S') B^B^ system that takes 
into account this IP constraint and a Gaussian constraint 
on the sum of the two B decay times {tcp + tta tr) t o be 
equal to 2 r^o with an uncertainty of ^/2 Tpo [a llO| . In 
order to ensure a well-determined vertex separation be- 
tween i?ioc and -Btag, we exclude events that have the 
error on At, determined from the vertex fit, a At > 2.5 ps 
and events with |At| > 20 ps. The mean uncertainty in 
Zcp, a convolution of the uncertainty in the interaction 
region position and the ztag resolution, is 75/im. The 
mean uncertainty on ztag is about 200/xm, which dom- 
inates the Az uncertainty. The resulting Az resolution 
is comparable to that in 5° J/ipK^ decays [sj. Sim- 
ulation studies and a 5° — + J/tpK'^ data control sample 
show that the procedure we use to determine the vertex 
for a Bcp decay provides an unbiased estimate of zcp . 

Most events have at least one candidate for which 
both tracks have at least one hit in the inner three SVT 
layers. We have verified on simulation and on data con- 
trol samples that the parameters of the signal At resolu- 
tion function for these Bcp signal decays are similar to 
those obtained from the -Bflav sample [3]. When at least 
one K'^ has tracks with hits in the outer two SVT layers 
but not in the inner three layers, the resolution is nearly 
two times worse and the At information is not used. 

We extract the event yields and CP parameters with an 
unbinned extended maximum-likelihood fit to the kine- 
matic, event shape, and At variables. For each of the 
sub-samples /c = 1,2 {Bcp{-^ )i^cp(oo)) use: 

Nt n 
i j 

where Vj is the probability density function (PDF) for 
the fit component. Nj is the event yield of each 
of the n components: Ns signal events, Nqq contin- 
uum qq events and, for Bcp{oo) only, Ngg BB back- 
ground events; Nt is the total number of events se- 
lected. For Bcp{-\ ) (-Bcp(oo)) candidates, the PDF 

Vj is given by the product of 'PjimEs)'Pj{AE)'PjiJ') 
iVj{m,,iss)Vj{mB)Vj{l2)) X V^iAt, aAt)s^, summed 
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over the tagging categories c. The product Hxlli is maxi- 
mized to determine the common CP asymmetry parame- 
ters S and C and the values of iVj , which are specific to 
each sub-sample. Along with S and C, the fit extracts 
e'^ and parameters describing the background. 

A fit to 857 Bcp(^ ) and 4992 -Bc'P(oo) candidates re- 
turns the event yields reported in Table|Tl Figure [1] shows 
the m-ES and AE (mmiss and uib) distributions for signal 
and background Bcp(^ ) (-Bcp(oo)) candidates. The ex- 
tracted CP parameters for the two separate sub-samples 
and the combined ones are shown in Table HI Using a 
Monte Carlo technique, in which we assume that the 
measured values for the CP parameters on the combined 
data sample are the true values, we find that the two sub- 
samples agree within 1.6a. The statistical significance of 
the CP violation is evaluated as a/2 • Aln(£i£2), where 
Aln(£i£2) is the change in the natural log of the com- 
bined likelihood for the no CP- violation hypothesis with 
respect to the maximum value. We estimate it to be 
2.9 standard deviations. Figure [2] shows distributions 



TABLE I: Event yields and CP asymmetry parameters ob- 
tained in the fit. The errors are statistical only. 




5.26 5.28 
mEs(GeV/c') 



0.1 
AE (GeV) 




5.3 5.4 

mg (GeV/c^) 



FIG. 1: Signal and background distributions of (a) rriES and 
(b) AE for S(7p(+_) candidates and of (c) mmiss and (d) ms 
for Bcp(oo) candidates. The signal and background distri- 
butions have been separated using the technique described 
in [ij . The curves represent the PDF projections. The back- 
ground distributions are shown in the insets. 



of At for and i?°-tagged events, and the asymmetry 
A(At) = (Nbo - iV^o) / (Nso + %o)- 

Systematic uncertainties on the CP parameters are 
given in Table HIl The systematic errors are evaluated us- 
ing large samples of simulated Bcp decays and the i?flav 
data sample. We perform fits to the simulated Bcp sig- 
nal with parameters obtained either from signal or i?flav 
events to account for possible differences in the At reso- 





Bcp(+~) 


BcP(QO) 


Combined 


Ns 


125 ± 13 


64± 12 




N,^ 


732 ± 28 


4942 ± 77 




Nbb 




-14 ±32 




S 
C 


-1 06 +°'^^ 

J^-^" -0.16 

— fl 08 +"•23 
U.UO _o 22 


0.24 ±0.52 
0.23 ±0.38 


-0.71 ± 0.24 
0.02 ± 0.21 




FIG. 2: Distributions of At for events weighted using the 
technique described in [ill] for Stag tagged as (a) or (b) _B°, 
and (c) the asymmetry ^(At). The points are the weighted 
data and the curves are PDF projections. 



lution function. We use the differences in the resolution 
function and tagging parameters extracted from these 
samples to vary the signal parameters. We account for 
possible biases due to the vertexing technique by com- 
paring fits to a large simulated sample of IP-constrained 
(neglecting the J/ip contribution to the vertex and us- 
ing the K2; trajectory only) and nominal J/ipK^ 
events. Several SVT misalignment scenarios are applied 
to the simulated Bcp events to estimate detector effects. 
We consider variations of 20/xm in the direction orthog- 
onal to the beam axis for the IP position and resolution 
and find they have a negligible impact. The systematic 
error due to correlations between the variables used in 
the fit is determined from a fit to a sample of randomly 
selected signal Monte Carlo (MC) events added to back- 
ground events generated from the background PDFs used 
in the fit. The values of the effective CP parameters 
for the BB background, which are fixed to zero in the 
nominal fit, are varied over the whole physically allowed 
range. The largest deviations in S and C resulting from 
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this variation are used as systematic uncertainties. The 
world-average values of Am^ and of the B'^ mean lifetime, 
T30, held fixed in the fit, are varied by their uncertain- 
ties 2] . We account for the possible interference between 
the suppressed b ucd and the favored b — > cud am- 
plitudes for some -Btag decays [l^]- Finally, we include 
a systematic uncertainty to account for imperfect knowl- 
edge of the PDFs used in the fit. Most of this uncertainty 
is due to MC statistics, the rest to differences between 
data control samples and MC simulation. 



TABLE II: Systematic uncertainties on S and C. 




a{S) 


a{C) 


Vertex reconstruction 


0.016 


0.003 


Resolution function 


0.005 


0.007 


Flavor tagging 


0.009 


0.015 


SVT alignment and IP position 


0.016 


0.008 


Fit correlation 


0.004 


0.025 


BB CP, Arud and Tgo 


0.008 


0.009 


Tag-side interference 


0.001 


0.011 


PDFs 


0.026 


0.031 


Total 


0.037 


0.046 



In summary, we measured the B° K°K°K° time- 
dependent CP asymmetries, S = -0.71 ± 0.24 ± 0.04 
and C = 0.02 ± 0.21 ± 0.05, where the first errors are sta- 
tistical and the second systematic. The statistical corre- 
lation between S and C is —14.1%. These results agree 
well with the SM expectation. This measurement, which 
is limited by the small statistics of the sample, constrains, 
but does not exclude contributions from physics beyond 
the SM, such as the low-energy supersymmetry [5]. These 
results supersede our previously pub lished CP asymme- 
try results for K",K°K° [13| and are consistent 
with the measurements performed by the Belle collabo- 
ration reported in [H . 
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